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ABSTRACT 

A methodology is developed to computationally simulate the uncertain 
behavior of composite structures. Uncertain behavior is the consequence of the 
random variation (scatter) of the primitive (independent random) variables at the 
constituent, ply, laminate and structural levels. This methodology is 
implemented in the IPACS (Integrated Probabilistic Assessment of Composite 
Structures) computer code. A fuselage-type composite structure is analyzed to 
demonstrate the code's capability. The probability distribution functions of 
structural responses are computed. Sensitivity of a given structural response 
to each primitive variable is also determined from the analyses. 


INIRDDUCnQN 


Composites are an important class of engineering material in aerospace 
design. Their outstanding mechanical properties, particularly the high strength 
to density ratio, are very attractive to the aerospace industry. They also 
posses excellent durability and corrosion-resistance. The mechanical properties 
of composite materials are derived from those of fiber and matrix constituents, 
which depend on a wide variety of variables . These variables include, for 
example, fiber and matrix material properties (constituent level) , fabrication 
variables such as fiber volume ratio, void volume ratio, ply orientation, ply 
thickness (ply level) . These primitive variables are statistical in nature, 
therefore, the mechanical properties of a composite material should be quantified 
probabilistically . As a consequence, the behavior of composite structures shews 
considerable scatter from its average value. 


Traditionally, a "safety factor" is used in designs to account for these 
unpredictable behavior. However, this approach results in either an 
overconservative design or an unsafe design. Because of the need in design for 
a probabilistic assessment capability for composite structures, a methodology was 
developed at the NASA Lewis Research Center. The methodology, incorporated into 
computer code IPACS, is based upon the identification of uncertain variables at 
every structural level. The uncertain variables are then filtered through an 
analyzer which combines the composite mechanics, structural mechanics and 
probability theory to determine the probability distribution functions of desired 
structural responses. 

These types of problems can normally be solved by the Monte Carlo 
simulation method. But, it is computationally expensive. In order to save 
computational time, the newly developed methodology uses special probability 
algorithms to provide an efficient and affordable way for the probabilistic 
assessment of uncertain composite structures. 

A fuselage-type composite structure is analyzed to demonstrate the code's 
capability. Buckling, natural frequency and static analyses are performed. The 
probability distributions for critical buckling load, first natural frequency, 
nodal displacement at free end, nodal strain and stress at fixed end are 
computed. The sensitivity of a given structural response to each primitive 
variable is also determined. 

PROBABILISTIC ASSESSMENT OF OCMPOSITE ST55JCIURES 

Probabilistic assessment of composite structures requires the 
identification of uncertainties within the model at the micro and macro levels. 
For example, at the constituent level, the uncertain primitive variables include 
the material properties of fiber and matrix. At the ply level, the uncertain 
primitive variables include fiber and void volume ratio, ply orientation and ply 
thickness. At the structural level, uncertainties in the loadings, temperature, 
geometry or boundary conditions may also be included. Each primitive variable 
is characterized by a characteristic probability distribution function and is 
input into an analyzer that combines composite mechanics, structural mechanics 
and probability theory. From the analysis , structural responses and material 
properties are obtained which are represented by probability distribution 
functions. The computational approach is illustrated schematically in figure 1. 
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The analysis modules which comprise the IPACS code are described briefly in the 
following section. 


OCMEOTER CODE IPACS 

IPACS is a computer code for probabilistic analysis of composite 
structures. It integrates several NASA in-house computer programs developed in 
recent years such as OOBSTRAN (Composite Blade STRuctural ANalyzer (ref. 1) , 
PICAN (Probabilistic Integrated Composite Analyzer) and PROCSAN (PRQbabilistic 
Composite Structural ANalyzer) . OOBSTRAN generates finite element models for 
composite structural analyses. PICAN determines the perturbed and probabilistic 
composite material properties at the ply and laminate level . PROCSAN determines 
the perturbed and probabilistic structural response at structural, laminate and 
ply levels. PICAN and PROCSAN use the same module, FPI (Fast Probability 
Integrator) (ref. 2, 3) , for probabilistic calculation. The architecture of THE 
IPACS software system is depicted in figure 2 . Details of PICAN, PROCSAN and FPI 
modules are explained in the following. 


PICAN MODULE 


PICAN has evolved from the deterministic composite mechanics code ICAN 
(Integrated Composite Analyzer) (ref. 4) . It performs a comprehensive multi- 
level perturbation and probabilistic analysis of composite material properties 
considering the uncertainties in the primitive variables at all composite scales. 
Perturbed material properties are computed at every level corresponding to each 
of the perturbed primitive variables. With this information, mathematical 
functions are developed to represent the relationships between the material 
properties and the primitive variables. The functions can be expressed by 

N ( 1 ) 

M p = a 0+ Y, a i x i 

i= 1 


if the relationship is linear, or by 
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( 2 ) 


Mp a o 


+ E a i 

7=1 


n 


+ L A 

i=l 


if the relationship is nonlinear, where Mp is any material property; n is the 
number of primitive variables; X^s are primitive variables; and a Q , a jf b j are 
constants. PICAN uses the probability distribution functions of the primitive 
variables, along with eq. (1) or (2) to calculate the probability distribution 
function for any material property . 

PROCSAN 

PROCSAN is a finite element module for the perturbation and probabilistic 
structural analysis considering all possible uncertainties including uncertain 
composite material properties. From laminate theory, the relationship between 
resultant forces and strains are shewn in equation (3) . 


A C 
C D 


( 3 ) 


where N are the resultant forces and M are the resultant moments; e are the in- 
plane strains and k are the curvatures . For a set of A, C and D matrices 
determined using perturbed primitive variables, the corresponding structural 
responses are computed by PROCSAN. Therefore, the structural responses and the 
primitive variables can be related by mathematical functions as expresses by 


R = a 0 


n 


+ E ^ x, 


( 4 ) 


if the relationship is linear, or by 
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( 5 ) 


n n 

r = *0 + £ a i x i + E b i x * 

i- 1 i= 1 

if the relationship is nonlinear, where R is structural responses. The 
probability distribution functions of structural responses will be similarly 
determined by probabilistic methods outline above. 


FPI MODULE 


The FPI module evolved from the probability algorithm FORM (First Order 
Second Moment Reliability Method) . After the perturbation study, a database is 
created. The relationship between the structural responses and the primitive 
variables can be determined numerically by FPI. Together with the probability 
distribution functions of primitive variables, FPI computes the cumulative 
distribution functions (cdf) of structural responses. 

IPACS SIMULATION 

The uncertain primitive variables at constituent, ply, laminate and 
structural levels are first identified before performing a probabilistic 
assessment of composite structures using IPACS. These variables are then 
selectively perturbed several times in order to create a database for the 
determination of the relationship between the desired structural response (or the 
desired material property) and the primitive variables. For a given set of 
perturbed primitive variable, micro mechanics theory is applied to determine the 
corresponding perturbed mechanical properties at the ply and laminate level. 
Laminate theory is then used to determine the perturbed A, C and D matrices in 
equation (3) . A finite element perturbation analysis is performed to determine 
the perturbed structural responses corresponding to the perturbed A, C and D 
matrices. This process is repeated until enough data is generated and the proper 
relationship between structural responses and primitive variables can be 
determined. 

With the known probabilistic distributions of the primitive variables and 
a relationship between the structural response and the primitive variables, fast 
probability integration is performed. For every discrete response value, a 
corresponding cumulative probability can be computed quickly by FPI. This 
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process is repeated until the cumulative distribution function can be 
appropriately represented. The probabilistic material properties at ply and 
laminate levels are also computed in the same way as for the structural 
responses. The output information from FPI for a given structural response 
includes its discrete cdf values, the coefficients to represent a cumulative 
distribution function, sensitivity factors of the primitive variables for 
probabilistic structural responses. 


NUMERICAL EXAMPLES AND DISCUSSIONS 


A cylindrical shell of dimensions shown in figure (3) is analyzed by IPACS. 
Four vertical stringers and three circumferential frames shown in figure 3 (a) are 
used to strengthen the shell-skin structure shown in figure 3 (b) . Three hundred 
and twelve nodes, including ninety six duplicate nodes, are used to model the 
structure for the finite element analysis. The graphite/epoxy composite 
structure is partially fixed at one end by translational and torsional springs 
and is free at the other end as shown in figure 4 (a) and 4 (b) . The average ply 
thickness for the shell skin is 0.005 inches with a [0,-45,45,90] s composite 
configuration. The average ply thickness for the stringers are 0.05 inches with 
10 zero degree plies. 

The primitive variables identified at constituent level are: (1) 17 fiber 
material properties, (2) 12 matrix material properties as listed in Table 1. At 
the ply level, the fabrication variables are: (1) fiber volume ratio, (2) void 
volume ratio, (3) ply orientation and (4) ply thickness. Their statistical data 
are listed in Table 2. At structural level, translational and torsional spring 
constants for boundary rigidity are modeled as uncertain parameters. 
Uncertainties are also assigned to those nodes on the free edge of the vertical 
stringers in their out-of-plane directions. Their statistics are shewn in Table 
3. Three probabilistic analyses are performed: (1) buckling analysis, (2) static 
analysis and (3) natural frequency analysis. Ninety three uncertain primitive 
variables are considered including sixty six variables from stringer and skin at 
the constituent and ply levels, three variables for the uncertain boundary 
conditions and twenty four out-of-plane nodal uncertainties in the vertical 
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stringers . 


BUCKLING ANALYSIS 

The probabilistic critical buckling load is computed. The cumulative 
distribution function (odf) of the buckling load is plotted in figure 5. The 
sensitivity factors of the six most important primitive variables for the 
critical buckling load at 0.005 cumulative probability are plotted in figure 6. 
It shows that the probabilistic buckling load is most sensitive to the ply 
thickness, followed by fiber volume ratio, fiber modulus of the skin and stringer 
ply thickness, respectively. The uncertainty in the torsional springs stiffness 
also has seme effect on the scatter of the probabilistic buckling load. 

FREQUENCY ANALYSIS 

Frequency analysis is also conducted. The mean value of the first natural 
frequency is 1.1 rad/sec. The cdf of this frequency is plotted in figure 7. The 
sensitivity factors of primitive variables for the natural frequency at 0.005 
cumulative probability is shewn is figure 8. Again, it is found that skin ply 
thickness is the most important variable with a sensitivity factor equal to 0.88. 

STATIC ANALYSIS 


Static analysis is performed using a load equal to the critical buckling 
load. The displacement at the free end in the longitudinal direction is 
analyzed. Its odf is plotted in figure 9. The sensitivity analysis shows that 
the fiber modulus, fiber volume ratio and ply thickness of both skin and stringer 
have approximately the same contribution to the scatter of the displacement. The 
sensitivity factors at 0.005 cumulative probability are plotted in figure 10. 
Also studied are the strain and stress of the node at the fixed end in the 
vertical stringer. The odf of the strain and its sensitivity factors are plotted 
in figures 11 and 12 respectively. The odf of stress and its sensitivity factors 
are plotted in figures 13 and 14 respectively. The sensitivity analysis again 
shews that the scatter of the strain and stress all are strongly dominated by the 
ply thickness, fiber volume ratio and fiber normal modulus of both skin and 
stringers. 
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OCNCUJDING REMARKS 


A methodology is developed to computationally simulate the uncertainties 
in the composite structural analysis. The methodology is implemented into the 
computer code IPACS which can efficiently and effectively simulate the uncertain 
structural responses. The code is capable of handling the uncertainties in the 
constituent, ply and laminate levels of composite structures. A large amount of 
useful information such as the probability distributions of desired structural 
responses and the sensitivity of the probabilistic structural responses to each 
independent primitive variable. 
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SYMBOLS 


E f11 : fiber modulus in longitudinal direction 
E f22 : fiber modulus in transverse direction 
G f12 : in-plane fiber shear modulus 
Gf 2 3 : out-of -plane fiber shear modulus 
v f12 : in-plane fiber Poisson's ratio 
v f23 : out-of-plane fiber Poisson's ratio 

a f11 : fiber thermal expansion coefficient in longitudinal direction 

a f22 : fiber thermal expansion coefficient in transverse direction 

p f : fiber mass density 

N f : number of fiber per end 

d f : filament equivalent diameter 
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fiber heat capacity 

fiber heat conductivity in longitudinal direction 

fiber heat conductivity in in-plane transverse direction 

fiber heat conductivity in out-of-plane transverse direction 

fiber tensile strength 

fiber compressive strength 

matrix elastic modulus 

matrix shear modulus 

matrix Poisson's ratio 

matrix thermal expansion coefficient 

matrix mass density 

matrix heat capacity 

matrix heat conductivity 

matrix tensile strength 

matrix compressive strength 

matrix shear strength 

matrix moisture coefficient 

matrix diffusivity 

fiber volume ratio 

void volume ratio 

ply misalignment 

ply thickness of skin 

ply thickness of stringer 

translational spring constant 

torsional spring constant 

out-of-plane nodal uncertainties on stringers at free edges 

standard deviation 

material properties 

structural response 

independent random variable 

strain 

curvature 

resultant in-plane force 

resultant moment 

number of primitive variables 
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Table l. Uncertain Material Properties at Constituent Level 
for both Skin and Stringers 



unit 

Distribution 

type 

Mean 

coefficient 
of variation 

E f11 

Msi 

Normal 

31.0 

0.05 

Ef22 

Msi 

Normal 

2.0 

0.05 

Gf12 

Msi 

Normal 

2.0 

0.05 

G f23 

Msi 

Normal 

1.0 

0.05 

V f 12 

— 

Normal 

0.2 

0.05 

V f23 


Normal 

0.25 

0.05 

“fll 

ppry°F 

Normal 

-0.55 

0.05 

a f22 

PPV F 

Normal 

5.6 

0.05 

Pf 

lb/ in 

Normal 

0.063 

0.05 

N f 

— 

constant 

10,000 

0.00 

d f 

in 

Normal 

0.0003 

0.05 

Cf 

BIU/in/°F 

Normal 

0.17 

0.05 

K fU 

BIU. in/hr/iny°F 

Normal 

580 

0.05 

K f22 

BIU. in/hr/ inV°F 

Normal 

58 

0.05 

K f33 

BIU. in/hr/ in / °F 

Normal 

58 

0.05 

S fT 

Ksi 

Weibull 

400 

0.05 

S fC 

Ksi 

Weibull 

400 

0.05 

E m 

Msi 

Normal 

0.5 

0.05 

G m 

Msi 

Normal 

0.185 

0.05 

V m 

— 

Normal 

0.35 

0.05 

a m 

ppm/°F 

Normal 

42.8 

0.05 

Pm 

lb/ in 5 

Normal 

0.0443 

0.05 

C m 

BIU/in/°F 

Normal 

0.25 

0.05 


BIU. iri/hr/in /°F 

Normal 

1.25 

0.05 

S mT 

Ksi 

Weibull 

15 

0.05 

S mc 

Ksi 

Weibull 

35 

0.05 

S mS 

Ksi 

Weibull 

13 

0.05 


in/in/l%moist 

Normal 

0.004 

0.05 

D m 

in /sec 

Normal 

0.002 

0.05 
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Table 2. Uncertain Fabrication variables at ply level 



unit 

Distribution 

type 

mean 

Coefficient 
of variation 

fvr 

— 

Normal 

0.60 

0.05 

wr 

— 

Normal 

0.02 

0.05 

e P 

degree 

Normal 

0.00 

0.9 (stdv) 

t psk 

in 

Normal 

0.005 

0.05 

tpst 

in 

Normal 

0.02 

0.05 




Table 3. Uncertain Variables at the 

structural 

level 



unit 

Distribution 

type 

mean 

coefficient 
of variation 


K CTR 

lb/in 

Normal 

30E+06 

0.30 


K cto 

lb- in/ rad. 

Normal 

12E+02 

0.30 


C z 

in 

Normal 

0.0 

0.01 (stdv) 
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COMPONENT 




COMPOSITE 

MICRO-MECHANICS 

THEORY 



CONSTITUENT 

PROPERTIES 


NONLINEAR 
MATERIAL BEHAVIOR 
MODEL 

<c= 


PLY RESPONSES 


COMPOSITE 

MICRO-MECHANICS 

THEORY 


DEGRADED 

MATERIAL PROPERTIES 


Figure 1 . Concept of Probabilistic Assessment of Composite Structures 



TABLES GRAPHICS 


Figure 2. Architecture of IPACS Software System 

(IPACS: Intergrated Probabilistic Assessment of Composite Structures) 
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(a) frames/stringers 


(b) skin 


Figure 3. Finite Element Model of a Composite Cylindrical Shell 



Figure 4. Boundary Conditions 



(b) torsional springs 

a Composite Cylindrical Shell 


Cumulative 

Distribution 

Function 



Figure 5. Cumulative Distribution function of first buckling load 
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Efl 1 of skin 


fvr of skin 

7 / 

skin ply thickness 


stringer ply thickness 


rotational spring Y-Y 

Z 

rotational spring Z-Z 


Figure 6. Sensitivity factors of six most important primitive 
variables for the first buckling load at probability level 0.005 
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Cumulative 

Distribution 

Function 



First Natural Frequency (rad/sec) 

Figure 7. Cumulative Distribution function of first natural frequency 


i 

0.8 

0.6 

Sensitivity 

Factors 

0.4 

0.2 
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Figure 8. Sensitivity factors of six most important primitive 
variables for the first natural frequence at probability level 0.005 
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0.25 
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Free End Displacement in X Direction 
(E-03 in) 


Figure 9. Cumulative Distribution function of free end displacement 
in X direction 


Sensitivity 

Factor 



Figure 10. Sensitivity factors of six most important primitive variables 
for the free end displacement in x direction at probability level 0.005 
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Cumulative 

Distribution 

Function 



(E-06) 


Figure 1 1 . Cumulative Distribution function of fixed end strain 
in X direction 
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fvr of skin 


& 

Ef 1 1 of stringer 

& 

fvr of stringer 
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Figure 12. Sensitivity factors of six most important primitive variables 
for the fixed end strain in x direction at probability level 0.005 
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0.75 - / 

Cumulative / 

Distribution 0.5 - / 

Function / 

0.25 - / 

0 I 1 ■ -" 1 ! 1 

-700 -630 -560 -490 -420 

Fixed End Stress in X Direction 
(psi) 

Figure 13. Cumulative Distribution function of fixed end stress 
in X direction 
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0.4 
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fvr of skin 


skin ply thickness 

A\ 

W 

Efl 1 of stringer 


fvr of stringer 

2 

stringer ply thickness 


Figure 14. Sensitivity factors of six most important primitive variables 
for the fixed end stress in x direction at probability level 0.005 
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